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Abstract 

We construct the cosmological model to explain the cosmological constant problem. We built the 
extension of the standard cosmological model ACDM by consideration of decaying vacuum energy 
represented by the running cosmological term. From the principles of quantum mechanics one can 
find that in the long term behavior survival probability of unstable states is a decreasing function of 
the cosmological time and has the inverse power-like form. This implies that cosmological constant 
Pvac = A(t) = Abare + ^ where Abare and a are constants. We investigate the dynamics of this model 
using dynamical system methods due to a link to the A{H) cosmologies. We have found the exact 
solution for the scale factor as well as the indicators of its variability like the deceleration parameter 
and the jerk. From the calculation of the jerk we obtain a simple test of the decaying vacuum in the 
FRW universe. Using astronomical data (SNIa, H{z), CMB, BAO) we have estimated the model 
parameters and compared this model with the ACDM model. Our statistical results indicate that 
the decaying vacuum model is a little worse than the ACDM model. But the decaying vacuum 
cosmological model explains the small value of the cosmological constant today. 
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I. INTRODUCTION 


The standard cosmological model (ACDM model) is a very good description of the evo¬ 
lution of the Universe in the current epoch. However it is only an effective description in the 
terms of dark matter and dark energy of unknown nature. The natural way of explanation 
why our Universe’s expansion is accelerating is to extend the CDM model with the constant 
cosmological term. On the other hand it seems to be natural to interpret the cosmological 
constant as the vacuum energy [1]. Although such an explanation brings the answer quali¬ 
tatively the calculations of the value of the cosmological constant term from the quantum 
field theory and its value required for solving the Universe accelerated expansion conun¬ 
drum give rise to discrepancy, unknown in the history of physics, of more than 100 orders 
of magnitude. Therefore the problem of cosmological constant is crucial for consistency of 
the standard cosmological model. 

In the literature different approaches have been taken to explain why the value of the 
cosmological constant is so small today. The most popular approach is description of the 
dark sector of the Universe in terms of the self-interacting scalar field with a potential [ 2 ]. 
Alternatively in this paper we develop the Alcaniz and Lima approach in which cosmology 
with the decay vacuum energy is interpreted in the framework of the A(f)CDM cosmology 
M. Alcaniz and Lima [3] pointed out that the critique of these models concentrates on 
lack of the theoretical model which can specify the phenomenological form of the relation 
describing how A depends on time. In this paper we are going to address the theoretical 
foundation of this approach. 

The discussion of false vacuum in the context of cosmology has been started by Coleman 
iEl- The decaying of false vacuum in the context of inflation has been investigated by 
Krauss and Dent [S] . The properties of unstable false vacuum can be investigated from the 
point of view of the quantum theory of unstable states p[in] and some universal behavior 
of survival probability of vacuum states can be discovered. 

From the quantum theory the amplitude A{t) and the decay law Viviit) of unstable state 
\M) are determined by energy density distribution function u:{E) for the system in this 
state. 

A{t) = [ u{E) exp (——Etj dt. (1) 

-'Emin V ^ / 

where uj{E) > 0 and uj{E) = 0 for U < 


2 


The last condition and the Paley-Wiener theorem give us |^(t)| > exp(—for 
t —)■ oo. In consequence the vacuum unstable states decay following the power-law relation 
after exponential for t S> T where T is the characteristic time of a regime switching. The 
asymptotically we obtain universal relation 

2 

± ^ ± ■ • • , for t > T. (2) 

At the late time survival probability assumes the form of the power-law relation. We 
demonstrate that then the instantaneous energy of the vacuum states has universal asymp¬ 
totics pvac = Pvac(t) = Pq™® + ^ Usiug this parameterization of the false vacuum 

will study the evolution of the Friedman-Robertson-Walker (FRW) universe hlled with dust 
matter and dark energy running in the cosmological time following the rule 

(y 

Pvac ■^(^) Ajj^j-e -|- — (3) 

where we use the natural system of units SnG = c = 1. 

The parametrization ([^ is only a leading term at the late time of evolution. In general 
from principles of quantum mechanics it assumes more general form in the form of power 
law series [9] 

Ht) = ^+ ^ + ^ + ^ + -" ( 4 ) 

The relation ([^ was proposed ad-hoc in the context of cosmological constant problem by 
Lopez and Nanopoulos HU. They showed that starting at the Planck epoch Api Mpj the 
Universe reaches today epoch with Aq ~ Another similar result was obtained by 

Lima et ah [T2] . 

The main aim of this paper is to explain the cosmological constant problem. We pro¬ 
pose the theoretical justification of choosing the varying A during the cosmic evolution 
parametrized as ([^. The form of this parametrization is motivated by quantum mechanics 
principles. The key point of this approach is the connection of decaying vacuum with mea¬ 
surement of a running jerk. The dimensionless parameter a is estimated from astronomical 
data. 
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II. THE EXACT SOLUTION FOR THE SCALE FACTOR IN THE FRW DECAY¬ 
ING VACUUM COSMOLOGY 

Throughout this paper we consider a flat homogeneous and isotropic universe in which a 
source of gravity is the perfect fluid with energy density p{t) and pressure p{t). We introduce 
the energy momentum tensor 

= + (5) 

where the first component describes dust and the second one the time-dependent A term. 
The energy-momentum conservation condition = 0 gives rise to the condition 

Pm T 3—Pjn Pvac (6) 

a 

where a(t) is the scale factor, pm and pvac are the energy densities of cold dark matter and 
vacuum, respectively, a dot means the differentiation with respect to the cosmological time t, 
Pvac = ^ and the equation of state for vacuum is Pvac = —Pvac where Pvac is pressure. 

The Einstein equations reduce to the acceleration equation 

3^ = -h„ + A(«) (7) 

(Jj Zd 

where A(f) is given by formula ([^. Equation ([^ can be interpreted formally within the 
framework of interacting cosmology 

Pm + 3-pm = r(f) (8) 

a 

where T = —Pvac = is the energy exchange term. This term is called a source term for 

the particle creation process [13]. If T(t) is positive then the vacuum energy density must 
be decaying. Because of parametrization ^ we have r(f) = 2at~^, which is positive for the 
positive a. 

For our aim it would be useful to rewrite the vacuum dark energy into the new 
parametrization. The idea of Alcaniz and Lima [3| is that the decay vacuum should modify 
the law of evolution of the cold dark matter (pm oc a~^) in such a way that matter dilutes 
more slowly like 

Pm ~ Pm, 0 ® (9) 

where e is a small positive constant in (0,1). 
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Thus, we substitute equation ([^ into the continuity condition ([^ and obtain 


Pvac PvaCjO T 


3 




— 3 - 1 -£ 


( 10 ) 


The constant Pvac,o is called “the ground state value of the vacuum” [1]. We interpret 
Pvac,o as Abare ~ the coustaut contribution to the vacuum energy parametrization. And we 
interpret the last term in equation (10) as a/t^ following the parametrization (|^. Hence 
after elementary calculations we obtain that 

a 3 — e 


-aH^ 


( 11 ) 


or 


A = A{H) = Abare + 3^H^ (12) 

where (3 = We assume that 0 < e < 1, and 

As a consequence we obtain the phenomenological quadratic parametrization of A through 
the Hubble parameter H [T3] . The same form of parameterization can be also obtained after 
replacing the cosmological time t by the Hubble scale tu = ^/H. 

Now it would be useful to rewrite the acceleration equation to the new form (where 
constant contribution Abare is denoted as A in further analysis) 


where 5 is constant defined as 


dt 2 


d=-{l-(3). 


(13) 


(14) 


It is convenient for our purposes to introduce a new time variable into equation (13). Let 


assume that there is the inverse function of t = t{a) for the cosmic evolution. Then 

dH dH da dH A 

dt da dt da 2 


(16) 


or 


da aH \ 2 


(16) 


where a prime denotes a derivative with respect to the scale factor a. Note that if /3 = 0 or 
6 = 3/2 we obtain a corresponding relation for the ACDM model. 


Now one can check that the solution of equation (16) is given in the form 


H^(a) ={H^-^\a 


-25 


A 


5 


(17) 








where H = Hq for a = oq = 1 and all quantities with a subscript “0” are evaluated at the 
present epoch. 

In the context of cosmography it would be useful to rewrite relation 0 to the new 
form in which appears dimensionless density parameters for the cold dark matter and the 
constant contribution to the dynamical dark energy, A, parameterized following formula ([^ 

= Jk- (18) 


^m,0 — 


Pm,0 

Wl' 


SHS' 


Then we obtain the relation H{z) for the extended ACDM model with decaying vacuum 
dark energy in the form 


H{z) 

Ho 


Q 


m,0 


n 


m,0 




(1 + m, 0+fiA,o) 




n 


m,0 


^A,C 


(19) 


where the redshift z = {1 + a)“^. If r2m o + ^A,o = 1 we obtain the corresponding relation for 
the ACDM model. Therefore, the effects of decaying vacuum are manifested by violation of 
the condition + ^a,o = 1- 

The parameter /3 in the parameterization of A{H) ( |I^ was constrained by using the 
observations of cosmic star formation and the age of the Universe as 3/3 < 0.3 [13]. In the 
next section we will search for more stringent upper limit for this parameter. 

Let us consider an exact solution for the function ([T^. The scale factor a satishes a hrst 
order differential equation 


a 2 _ tt2 

a — Uq 


n 


m,0 


^m ,0 “ 1 “ ^Ao 


j —3(Qm,0+^A,o)H“2 I 




^m,0 “1 “ ^Aq 


( 20 ) 


Equation (20) can be integrated after a simple substitution 


a ^ x: = a3(nin,o+OA,o) 


( 21 ) 


and the time parametrization 


t^T: T = -|J?oI\/11a,o<' 


( 22 ) 


After solving the equation for x{t) and returning to original variables we obtain the following 
expression for the scale factor as a function of the cosmological time 


a{t) = 


Dm,o\ 8(nin,0+flA,o) 

f^A,0 


sinh -^Da,o 


3(llm,0+flA,o) 


(23) 
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Formula (23) gives us the possibility to calculate the cosmological scalars which control the 
variability of the cosmic evolution 



a dt 



(24) 


These scalars are called the Hubble parameter, deceleration parameter, jerk, snap and 
crackle, respectively. They are constrained by some algebraic relation [T5] . 

Let us calculate the cosmic jerk related with a third order time derivative of the scale 
factor 



(25) 


2 

3 


where A = 

If A = I we have the case of the ACDM model for which Q = 1 and if A 7 ^ 



Hm o + ^A,o then Q = Q{t). Therefore if we hnd that the jerk is different from one it 
will be experimentum crucis for excluding the ACDM model. Then as the alternative the 


cosmological model with decaying vacuum should be treated seriously. In this model the 


jerk is time dependent and now different from one as the function of sum of Dm,o and f 2 A,o 
which is different from one. 

III. DYNAMIC OF THE MODEL ON THE PHASE PLANE 

For deeper analysis of the dynamics of the model let us formulate this model as a dy¬ 
namical system. 


h,„ + i(A + A(i/)) 

0 3 

p„ = -3Fp„ + j(A + k{H)) 


H = 


(26) 


(27) 


7 








Dynamical system (26)-(27) assumes the form of two-dimensional autonomous dynamical 
system. In the special case of A.{H) = aH^ and 






<y „ /a ,^3 2 


Pm = ( 2 - 3j if Pm + 2a - Ij H - -aH. 


(28) 

(29) 


It could be useful to rewrite the above system in two state variables and new time variable 


x = H^, y = pm, T = \na. 


(30) 


This system possesses the critical point 


A 


Xo = 


Vo = 0. 


a-3’ 

Now we translate the critical point at the origin of the coordinate system as follow 

X = x-Xo Y = y-yo. 

Finally the dynamical system possesses the form of linear system 

A-V) ^ - ir(.) 

Y'ir) = Y±Jz3^X{r)+^-^Y(r) 


(31) 


(32) 


(33) 

(34) 


The phase portrait for system (33)-(34) with the positive cosmological constant is pre¬ 


sented in Fig. The system possesses the critical point which is stable node. It is rep¬ 
resenting the de Sitter solution with the cosmological constant A. The physical region is 
determined by condition x = < 0 and p > 0. Therefore the physical trajectories should 

lie in the hrst quarter of (x, p)-coordinate systems. Note that if A is negative the unique 
critical point is situated in the non-physical region. 

The solution of the system in the parametric form is 

1 


X(r) = 
Y{t) = 


3(a- 1) 
1 


-1- (a - 3)6-^") Ci + Ca] 


(35) 


[0 


,(«- 3 )t _ - 2 t 


e-2") 2a(a - 3)^1 + {2ae-^^ + {a - 3)e(“-^)") C 2 ] (36) 


3(a- 1) 

Note that there are four special solution lying along two the straight lines 


Y{t) = f)A(r) 


(37) 












0 



FIG. 

A{H) 


1. The phase portrait of system (33)-(34) with the positive cosmological constant A and 
= aH^. The critical point at the origin represents the de Sitter solution. 


where the parameter /3 is 


or 


(3 = — a 

13 = 2a. 


(38) 

(39) 


The line a; = 0 separates the cosmological models with a singularity from singularity-free 
models. The latter are starting from the boundary of non-physical region and going toward 
the de Sitter global attractor. The former are starting from singularities {H = oo, = oo) 
of the Einstein-de Sitter or Milne type. 

For full analysis of the system dynamics it is necessary to know how the system trajectories 
behave at the inhnity. For this aim we will construct the global phase portraits. In the case 
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when right-hand sides of a dynamical system are given in a polynomial form of order n the 
projective coordinates can be introduced. For a two-dimensional dynamical system with 
state variables (x, y) two projective maps 

I: z = -, u=- (40) 

X X 

II; w = -, V = —. (41) 

y y 

cover the dynamics. Of course when m 7^ 0 and n 7^ 0 these two projective coordinates {z, u) 
and {v,w) are equivalent. The inhnitely distant points of the plane {x,y) in the projective 
coordinates correspond to a circle : x^ + y‘^ = 00 covered by two straight lines: z = 0, 
—00 < u < -|-oo and tc = 0 , —00 < v < -|-cxd. The original dynamical system in state 
variable {x,y): x = P{x,y), y = Q{x,y) after transforming to the projective coordinates 
and reparametrization of the time variable r —)■ Ti : dxi = x'^~^dT assumes the form 


where 


z = —zP*{z, u), 
u = Q* {z, u) — uP* {z, u) 


p-{z, u) = ( i, - j , Q’(z,u) = 

ZZ \Z Z 


(42) 

(43) 

(44) 


are polynomials in the projective variable {z,u). 

The analogous procedure can be performed in the map (n, w) 


Dynamical system (33)-(34) in the projective coordinates z{t) = 7 ^, u{t) = has 
the following form 


z')’’) = ]^ u { t ) z { t ) - ?(a - 3 ) z ( t ) 

'“'(t) = “ ^(“ + 3)“('r) “ 


(45) 

(46) 


The exact solution of the system (45)-(46) assumes the following form 


Z T = 


U T = 


3(a-l) 


[2ae(“-3)^ + (a- 3)e-2^] Ci + (e-^^ - e("-3)^) C 2 
(e(“-3)^ _ e-2^) 2a(a - 3)Ci + (2ae-^^ + (a - 3)e("-3)^) C 2 


(47) 

(48) 


[2ae(“-3)^ (a - 3)e-2^] Ci + (e- 2 ^ - e("-3)^) C 2 

The critical points are z = 0, u = 3 — a and z = 0, u = 2a. The former is an unstable node 
and the latter is a saddle. 
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Dynamical system (33)-(34) in the projective coordinates w{t) = ~ 


the following form 


2 1 

w\t) = -{a — 3)Q;tc(r)n(r) — -(a — 9)w{t) 
o o 

= j(a - 3)av^(T) + ^(a + 3)v(t) - t. 


(49) 

(50) 


The exact solution of the system (49)-(50) assumes the following form 

3(a- 1) 


w(t) = 
v(t) = 


(^e(a 3 )t _ g 2 tj 2a(a — 3)Ci + [206 + (a — 3)e(" C 2 

(2ae(“-3)^ + (a - 3)e-2^) Ci + (e-^^ - C 2 

(e(““3)T _ e- 2 T^ 2a{a — 3)C'i + [2Q;e“^’' + {a — 3)e("“^)’'] C 2 


(51) 

(52) 


The critical points are w = 0, v = and w = 0, v = ^. The former is an unstable node 
and the latter is a saddle. 

The phase portrait of the system under consideration are presented in Fig. an di From 
these phase portraits one can observe how trajectories behave near the initial singularity. In 
both cases two critical points are situated on the lines z = 0 and n = 0. The phase portrait 
is structurally stable which means that the system is resistant with respect to small changes 
of its right-hand sides. 

The above analysis can be continued for cosmological models with a general form A = 
A{H). Recently Perico et ah have studied the dynamics and exact solution if A{H) is given 
in power series with respect to H [16]. The physical motivation for this studies can derived 
from quantum mechanics principles im. 


IV. DATA 

To estimate the parameters of the both models we used the modified for our purposes 
CosmoMC code [181 [T9| with the implemented nested sampling algorithm multinest [201 

m\. 

We used the observational data of 580 supernovae type la (the Union2.1 compilation 
|22j). 31 observational data points of Hubble function [23], the measurements of BAO (barion 
acoustic oscillations) from the Sloan Digital Sky Survey (SDSS-III) combined with the 2dF 
Galaxy Redshift Survey (2dFGRS) [241127] . the 6dF Galaxy Survey (6dFGS) [28l |29], the 
WiggleZ Dark Energy Survey [301132] . We also used information coming from determinations 
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FIG. 2. The phase portrait of the system in the projective coordinates (z, u) dedicated for the 
analysis of the behavior of trajectories at the infinity —)• oo. There are two critical points: 

a saddle point representing the Milne universe oc a~^ and an unstable node representing the 
Einstein-de Sitter-like universe. It is assumed the parameter e is rather small. The shadowed 
region is representing the non-physical domain forbidden for trajectories where < 0. 


of the Hubble function using the Alcock-Paczyhski test [33l l8l] . This test is very restrictive 
in the context of modified gravity models. 

We use the following likelihood functions for aforementioned data. 

First, the likelihood function for the supernovae type la data 


Lsn oc exp 




(53) 


where is the covariance matrix with the systematic errors, — M is the distance 
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FIG. 3. The phase portrait of the system in the projective coordinates {v, w) dedicated for the 
analysis of the behavior of trajectories at the infinity —)• oo. There are two critical points: 

a saddle point representing the Milne universe oc a~^ and an unstable node representing the 
Einstein-de Sitter-like universe. It is assumed the parameter e is rather small. The shadowed 
region is representing the non-physical domain forbidden for trajectories where < 0. 


modulus, = 5 logio ^Li + M = 'b log^o + 2h, M = -5 log^o -^o + 25 and Du = Hodu, 

where dn is the luminosity distance which is given by dn = (1 + Zi)c fj" (with the 
assumption fc = 0). 

Second, the likelihood function for H{z) data 


oc exp 


E 


{H^'^izi) - 

2^1 


(54) 


where H^^{zi) denotes the theoretically estimated Hubble function, is observational 
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data. 


Third, the likelihood function for the BAO data 


Lbao oc exp 






(55) 


where Cij is the covariance matrix with the systematic errors, = rs{zd) (1 + Zi)‘^D\{zi)^j^^ 

rs{zd) is the sound horizon at the drag epoch and Da is the angular diameter distance. 

Fourth, the likelihood function for the information coming from the Alcock-Paczyhski 
test 

{AP^'^{zi) - 


Lap cc exp 


-E 


2<t? 


(56) 


where ^ 


Fifth, the likelihood function for the CMB shift parameter R 

1 - R^'^^y 


Lcmb oc exp 


a 


A 


(57) 


where R^'^ = -|_ z^DAiz^, DAy) is the angular diameter distance to the last 

scattering surface, R°^^ = 1.7477 and = 48976.33 [36] . 

And finally, the total likelihood function Ltot 


Ltot — TsNTn.TRAnT 


SN CMB J^AP- 


(58) 


V. STATISTICAL ANALYSIS 


A. Estimation of model parameters 


Let us assume that we have N pairs of measurements {yt, Xj) and that we want to find the 
relation between the y and x variables. Suppose that we can postulate k possible relations 
y = fi{x, 6), where 6 is the vector of unknown model parameters and i = 1,... ,k. With 
the assumption that our observations come with uncorrelated Gaussian errors with a mean 
/ij = 0 and a standard deviation ai the goodness of fit for the theoretical model is measured 
by the quantity given by 


N 


^' = E 


2=1 


{fi{xi,9) - yif 


= —2 InL, 


(59) 
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TABLE I. Mean of marginalized posterior PDF with 68% confidence level for the parameters 
of the models. In the brackets are shown parameter’s values of joined posterior probabilities. 
Estimations were made using the Union2.1, h{z), BAO, determinations of Hubble function using 
Alcock-Paczyhski test and CMB R data sets. 



Union2.1 data 

Union2.1, h{z), BAO, A-P 

Union2.1, h{z), BAO, A-P, CMB 

testing model 

^m,0 £ [0» 1] 

^^A,o e [0,1] 

hioo G [0.60, 0.80] 

0.3032(0.2916) 
0.6719l](;j((:53 (0.7069) 
0.6998l°;J°“ (0.7395) 

0.2880to;ol?%0.2814) 

0.6968to.oM%0.7059) 

0.6917t°;°°®®(0.6955) 

0.2798lQ;Qm%0.2818) 

0.7177l“;°“®2(o.7i59) 

0.6927lo;oo7| (0.6897) 

ACDM model 

llm,0 £ [0, 1] 

^100 £ [0.60, 0.80] 

0.2953l°;“°|® (0.2954) 
0.7004l°;'j'®®®(0.7789) 

0.2777tg;K^?^ (0.2795) 
0.6932t°°°®°(0.6917) 

0.2908l!(;g!(^| (0.2906) 

0.68601q;pq4%0.6868) 


where L is the likelihood function. For the particular family of models fi the best one 
minimize the quantity, which we denote fi{x,6). The best model from our set of k 
models fi{x,9),... ,fk{x,9) could be the one with the smallest value of quantity. But 
this method could give us misleading results. Generally speaking for more complex model 
the value of x^ is smaller, thus the most complex one will be choose as the best from our 
set under consideration. 


B. The likelihood ratio test 


We perform the likelihood ratio test at the beginning of our statistical analysis. In this 
test one of the models (null model) is nested in a second model (alternative model) by hxing 
one of the second model parameters. In our case the null model is the ACDM model, the 
alternative model is the decaying vacuum ACDM model, and the parameter in question is 
f^vac- 


T 1 

Hi : Dm + Da 7^ 1. 


The statistic is given by 


A = 2 In 


( HHAD) 


= 2 


A vac 


Xacdm ^ 

2 ) 


(60) 


where L{Hi\D) is the likelihood of the decaying vacuum model, L{Hq\D) is the likelihood 
of the ACDM model. The statistic A has the distribution with df = ni — Uq = 2 degree 
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TABLE II. The results of the likelihood ratio test for the ACDM model (null model) and the 
decaying vacuum model (alternative model). The values of Xtest) TacdM’ statistic A and 
corresponding p-values (d/ = 4 — 2 = 2). Estimations were made using the Union2.1, h{z), BAO, 
determinations of the Hubble function using Alcock-Paczyhski test, and CMB R data sets. 


data sets 

X?ac/2 

Aacdm/^ 

A 

p-value 

Union2.1 

272.5544 

272.5552 

0.0016 

0.9992 

Union2.1, h{z), BAO, AP 

282.2388 

282.2555 

0.0334 

0.9834 

Union2.1, h{z), BAO, AP, CMB 

282.2489 

282.4913 

0.4848 

0.7847 


of freedom where rii is number of the parameters of the alternative model, no is number of 
the parameters of the null model. The results are presented in Table |IT1 The three different 
sets of data have been used (the Union2.1 only, all data without CMB and all data). In all 
three cases the p-values are greater than the significance level a = 0.05, that why the null 
hypothesis cannot be rejected. In other words we cannot reject the hypothesis that there is 
no decaying vacuum. 


C. Bayesian analysis of the models 

For the further comparison of the decaying vacuum model with the ACDM model we 
use the Bayesian statistical methods. These methods of models selection are widely used 
for cosmological model comparison We should be aware that conclusions based on 

such analyses depend on the data at hand. New data obtained in future may change our 
conclusions. 

First we calculate the Akaike information criterion (AIC) |16] which is an approximation 
to the Kulback-Leibler information which measures information lost when an unknown (hy¬ 
pothetical) true model is approximated by a tested model. The AIC quantity is lowest for 
a best model and is given by 

AIC =-2 In £ + 2(i, (61) 

where C is the maximum of the likelihood function and d is the number of model parameters. 
It is convenient to evaluate the differences between the AIC quantities computed for the rest 
of models from the set and the AIC for the best one. Those differences (Aaic) are interpreted 
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as a ‘strength of evidence’ for a model considered with respect to the best one (i.e. with a 
lowest value of AIC). The models with 0 < Aaic < 2 have substantial support (evidence), 
those where 4 < Aaic A 7 have considerably less support, while models having Aaic > 10 
have essentially no support with respect to the best model. 

The complexity of the model is interpreted here as the number of its free parameters that 
can be adjusted to £t the model to the observations. If models fit the data equally well 
according to the Akaike rule the best one is with the smallest number of model parameters 
(the simplest one in such an approach). 

In the Bayesian framework the best model (from a model set under consideration) is that 
which has the largest value of probability in the light of data (so called posterior probability) 

m 


P{M,\D) = 


P(D|M,)P(M,) 

P(P) 


(62) 


where P{Mi) is a prior probability for the model Mi (it is assumed that there is no evidence 
to favor one model over another, so values of priors for all models under consideration are 
equal), D denotes data, P{D) is the normalization constant 

k 


P(P) = 5^P(P|M,)P(M,). (63) 

i=l 

And P{D\Mi) is the marginal likelihood, also called evidence 


P{D\M,) = j P{D\e,Mi)P{e\Mi) de = P„ (64) 

where P{D\9, Mi) is likelihood under model i, P{6\Mi) is prior probability for 6 under model 

i. 


It is convenient to evaluate the posterior ratio for models under consideration which in 
the case with flat prior for models is reduced to the evidence ratio called the Bayes factor 

P{D\Mi) 


P(D\Mj 


(65) 


The interpretation of twice the natural logarithm of the Bayes factor is as follow: 0 < 
21nPjj < 2 as a weak evidence, 2 < 21nPjj < 6 as a positive evidence, 6 < 21nPjj < 10 as 
a strong evidence and 2 In Bij > 10 as a very strong evidence against model j comparing to 
model i. Let us simplify the problem to illustrate how this principle works here [39l 148] . 

Assume that P{9\D,M) is the non normalized posterior probability for the vector 9 of 
model parameters. In this notation P = j P{9\D,M)d9. Suppose that posterior has a strong 
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peak in the maximum: 6*mod- It is reasonable to approximate the logarithm of the posterior 
by its Taylor expansion in the neighborhood of ^mod so we hnished with the expression 


Pie\D,M) = PieMOD\D,M)x 

X exp [—{6 — Omod)'^C ^{6 — ^mod)] , (66) 


■<-11 




In P{e\D,M) 


aeidOi 


where [C 
distribution with the mean 6'mod 
form 


. The posterior is approximated by the Gaussian 

d=SyiOT) 

and the covariance matrix C. The evidence then has a 


E = P{eMOJ^\D,M) X 

X y exp [-(0 - 0mod)^G"^(i9 - 0mod)] dO. (67) 

Because the posterior has a strong peak near the maximum, the most contribution to the 
integral comes from the neighborhood close to ^mod- Contribution from the other region of 9 
can be ignored, so we can expand the limit of the integral to whole With this assumption 
one can obtain = (27r)i\/det GP(i9mod|-D, M) = (271)^^/det C'P(i7|0MOD, dT)P(0MOD|dT). 
Suppose that the likelihood function has a sharp peak in 9 and the prior for 9 is nearly flat 
in the neighborhood of 0. In this case 9 = 9mod and the expression for the evidence 
takes the form E = £(27r)^\/det CP{9\M). The quantity {27i)^Vdet CP{9\M) is called the 
Occam factor (OF). When we consider the case with one model parameter with a flat prior 
P{9\M) = ^ the Occam factor OF= ^ which can be interpreted as the ratio of the volume 
occupied by the posterior to the volume occupied by prior in the parameter space. The more 
parameter space wasted by the prior the smaller value of the evidence. It is worth noting 
that the evidence does not penalize parameters which are unconstrained by the data |41j . 

As the evidence is hard to evaluate an approximation to this quantity was proposed by 
Schwarz |19] so called Bayesian information criterion (BIG) and is given by 

BIG =-2 In/: + 2d In AT, (68) 


where N is the number of the data points. The best model from a set under consideration 
is this which minimizes the BIG quantity. One can notice the similarity between the AIG 
and BIG quantities though they come from different approaches to model selection problem. 
The dissimilarity is seen in the so called penalty term: ad, which penalize more complex 
models (complexity is identihed here as the number of free model parameters). One can 
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TABLE III. Values of AIC, BIC and Bayes factor, AAIC (with respect to the ACDM model). 
Estimations were made using the Union2.1, h{z), BAO, determinations of Hubble function using 


Alcock-Paczyhski test and CMB R data sets. 



XV2 

Evidence In Ei 

AIC 

AICj - AlCi 

BIC 

2 In Bij 

tested model 

Union2.1 

272.5544 

— 275 1Q31+0-0847 

551.1088 

1.9984 

564.1979 

6.3614 

Union2.1, h{z), BAO, A-P 

282.2388 

-289.0193to°|8? 

570.4776 

1.9666 

583.7716 

6.3979 

Union2.1, h{z), BAO, A-P, CMB 

282.2489 

-292.4944l[J;5^g^ 

570.4978 

1.5152 

583.7966 

5.9481 

ACDM model 

Union2.1 

272.5552 

-274.7896i°;“i° 

549.1104 

— 

557.8365 

— 

Union2.1, h{z), BAO, A-P 

282.2555 

-287.1432t0;“52 

568.5110 

— 

577.3736 

— 

Union2.1, h{z), BAO, A-P, CMB 

282.4913 

-287.9753l°i"| 

568.9826 

— 

577.8485 

— 


evaluated the factor by which the additional parameter must improve the goodness of £t to 
be included in the model. This factor must be greater than a so equal to 2 in the AIC case 
and equal to InN in the BIC case. Notice that the latter depends on the number of the 
data points. 

It can be shown that there is the simple relation between the BIC and the Bayes factor 


2\nBij = -(BICi-BICj). 


(69) 


The quantity Bij is the Bayes factor for the hypothesis (model) i against the hypothesis 
(model) j. We categorize this evidence against the model j taking the following ranking. 
The evidence against the model j is not worth than bare mention when twice the natural 
logarithm of the Bayes factor (or minus the difference between BICs) is 0 < 2 In B^j < 2, is 
positive when 2 < 2 In Bij < 6, is strong when 6 < 2\nBij < 10 and is very strong when 
2\nBij > 10. 


To obtain the values of AIC and BIC quantities we perform the = ~2 In L minimization 
procedure after marginalization over the Hq parameter in the range [60,80]. They are 


presented in Table III 


VI. CONCLUSION 


The main goal of the paper was to hnd the theoretical explanation of the cosmological 
constant problem. To solve the disrepancy of the value of the cosmological constant inter¬ 
preted as vacuum energy and value of cosmological constant required for the explanation of 
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FIG. 4. Posterior constraints for testing model. Solid lines denote 68% and 95% credible intervals of 
fully marginalized probabilities, the colors illustrate mean likelihood of the sample. Top left: esti¬ 
mations with the Union2.1 data only. Top right: estimations made using the Union2.1, h{z), BAG, 
determinations of Hubble function using Alcock-Paczyhski test data sets. Bottom: estimations 
made using the Union2.1, h(z), BAG, determinations of Hubble function using Alcock-Paczyhski 
test and CMB R data sets. 

acceleration of the current Universe we construct the cosmological model with the decay¬ 
ing vacuum. The parametrization of running cosmological constant is taken directly from 
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FIG. 5. Posterior constraints for ACDM model. Solid lines denote 68% and 95% credible intervals 
of fully marginalized probabilities, the colors illustrate mean likelihood of the sample. Top left: esti¬ 
mations with the Union2.1 data only. Top right: estimations made using the Union2.1, h{z), BAG, 
determinations of Hubble function using Alcock-Paczyhski test data sets. Bottom: estimations 
made using the Union2.1, h{z), BAO, determinations of Hubble function using Alcock-Paczyhski 
test and CMB R data sets. 


quantum mechanics principles which distinguish the power-law relation for decaying part of 
vacuum dark energy. 

Let us compare both the value of vacuum dark energy at the Planck epoch and today. 
Then we obtain 

Pa(^ - ^Pl) _ ^ .ypN 

pK{t ^ to) 

Expressing Hq in the Planck unit {Hq = 1.24 ■ we obtain 


Phit ~ tpi) 
pK{t ~ to) 


~ a - 4.5 ■ 10^2^ 


(71) 


From the best ht value 1 — Gm,o ~ = 0.0025 we found that 


a = 3(1 - - ^^a) = 0.0074 ~ 10“^ 


(72) 
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and finally 


pkii — ^pi) ^ ini20 
pA{t ~ to) ~ 


(73) 


On the gronnd of onr statistical analysis we can conclnde that the cosmological model 
with decaying vacnnm shonld be rejected and the ACDM model can keep a crown of the 
standard cosmological model. We showed that the cosmological model with decaying vacnnm 
is worst-htted than the ACDM model, bnt we think that it shonld not be considered yet 
as a blind alley. The strength of this model is its explanation of the cosmological constant 
connndrnm. That why it shonld be considered serionsly as a best alternative for the ACDM 
model from the physical point of view. 

We share George Efstathion’s opinion [301 - 1^^ that withont sonnd theoretical basis there 
is no reason for considering the dynamical dark energy. This paper is the attempt to base 
the cosmological model with varying A on the qnantnm mechanics principles. We snbmit 
this model for consideration althongh it is worse than standard cosmological model in the 
statistical terms. Bnt it solves the cosmological constant problem. The ACDM model has 
only the status of the satisfactory effective theory. Our model has a deeper explanatory 
value and we think it indicates the direction of the future research in cosmology. 
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